Genome-wide co-expression analysis in multiple tissues

Expression quantitative trait loci (eQTLs) represent genetic control points of gene expression, and can be categorized as cis- and trans-acting, reflecting local and distant regulation of gene expression respectively. Although there is evidence of co-regulation within clusters of trans-eQTLs, the extent of co-expression patterns and their relationship with the genotypes at eQTLs are not fully understood. We have mapped thousands of cis- and trans-eQTLs in four tissues (fat, kidney, adrenal and left ventricle) in a large panel of rat recombinant inbred (RI) strains. Here we investigate the genome-wide correlation structure in expression levels of eQTL transcripts and underlying genotypes to elucidate the nature of co-regulation within cis- and trans-eQTL datasets. Across the four tissues, we consistently found statistically significant correlations of cis-regulated gene expression to be rare (<0.9% of all pairs tested). Most (>80%) of the observed significant correlations of cis-regulated gene expression are explained by correlation of the underlying genotypes. In comparison, co-expression of trans-regulated gene expression is more common, with significant correlation ranging from 2.9%-14.9% of all pairs of trans-eQTL transcripts. We observed a total of 81 trans-eQTL clusters (hot-spots), defined as consisting of > or =10 eQTLs linked to a common region, with very high levels of correlation between trans-regulated transcripts (77.2-90.2%). Moreover, functional analysis of large trans-eQTL clusters (> or =30 eQTLs) revealed significant functional enrichment among genes comprising 80% of the large clusters. The results of this genome-wide co-expression study show the effects of the eQTL genotypes on the observed patterns of correlation, and suggest that functional relatedness between genes underlying trans-eQTLs is reflected in the degree of co-expression observed in trans-eQTL clusters. Our results demonstrate the power of an integrative, systematic approach to the analysis of a large gene expression dataset to uncover underlying structure, and inform future eQTL studies

Health Quality Control of Laboratory Anjmals 0f the Czechoslovak Academy of Sciences (CSAV)

No abstract availabl

Cap analysis of gene expression reveals alternative promoter usage in a rat model of hypertension

The role of alternative promoter usage in tissue-specific gene expression has been well established, however, its role in complex diseases is poorly understood. We performed cap analysis of gene expression (CAGE) sequencing from the left ventricle (LV) of a rat model of hypertension, the spontaneously hypertensive rat (SHR), and a normotensive strain, Brown Norway (BN) to understand the role of alternative promoter usage in complex disease. We identified 26,560 CAGE-defined transcription start sites (TSS) in the rat LV, including 1,970 novel cardiac TSSs. We identified 28 genes with alternative promoter usage between SHR and BN, which could lead to protein isoforms differing at the amino terminus between two strains and 475 promoter switching events altering the length of the 5’ UTR. We found that the shift in Insr promoter usage was significantly associated with insulin levels and blood pressure within a panel of HXB/BXH recombinant inbred rat strains, suggesting that hyperinsulinemia due to insulin resistance might lead to hypertension in SHR. Our study provides a preliminary evidence of alternative promoter usage in complex diseases

New insights into the genetic control of gene expression using a Bayesian multi-tissue approach.

The majority of expression quantitative trait locus (eQTL) studies have been carried out in single tissues or cell types, using methods that ignore information shared across tissues. Although global analysis of RNA expression in multiple tissues is now feasible, few integrated statistical frameworks for joint analysis of gene expression across tissues combined with simultaneous analysis of multiple genetic variants have been developed to date. Here, we propose Sparse Bayesian Regression models for mapping eQTLs within individual tissues and simultaneously across tissues. Testing these on a set of 2,000 genes in four tissues, we demonstrate that our methods are more powerful than traditional approaches in revealing the true complexity of the eQTL landscape at the systems-level. Highlighting the power of our method, we identified a two-eQTL model (cis/trans) for the Hopx gene that was experimentally validated and was not detected by conventional approaches. We showed common genetic regulation of gene expression across four tissues for ∼27% of transcripts, providing >5 fold increase in eQTLs detection when compared with single tissue analyses at 5% FDR level. These findings provide a new opportunity to uncover complex genetic regulatory mechanisms controlling global gene expression while the generality of our modelling approach makes it adaptable to other model systems and humans, with broad application to analysis of multiple intermediate and whole-body phenotypes

Towards precision medicine for hypertension: a review of genomic, epigenomic, and microbiomic effects on blood pressure in experimental rat models and humans

Compelling evidence for the inherited nature of essential hypertension has led to extensive research in rats and humans. Rats have served as the primary model for research on the genetics of hypertension resulting in identification of genomic regions that are causally associated with hypertension. In more recent times, genome-wide studies in humans have also begun to improve our understanding of the inheritance of polygenic forms of hypertension. Based on the chronological progression of research into the genetics of hypertension as the "structural backbone," this review catalogs and discusses the rat and human genetic elements mapped and implicated in blood pressure regulation. Furthermore, the knowledge gained from these genetic studies that provide evidence to suggest that much of the genetic influence on hypertension residing within noncoding elements of our DNA and operating through pervasive epistasis or gene-gene interactions is highlighted. Lastly, perspectives on current thinking that the more complex "triad" of the genome, epigenome, and the microbiome operating to influence the inheritance of hypertension, is documented. Overall, the collective knowledge gained from rats and humans is disappointing in the sense that major hypertension-causing genes as targets for clinical management of essential hypertension may not be a clinical reality. On the other hand, the realization that the polygenic nature of hypertension prevents any single locus from being a relevant clinical target for all humans directs future studies on the genetics of hypertension towards an individualized genomic approach

Targeting of the Plzf gene in the rat by transcription activator-like effector nuclease results in caudal regression syndrome in spontaneously hypertensive rats

Recently, it has been found that spontaneous mutation Lx (polydactyly-luxate syndrome) in the rat is determined by deletion of a conserved intronic sequence of the Plzf (Promyelocytic leukemia zinc finger protein) gene. In addition, Plzf is a prominent candidate gene for quantitative trait loci (QTLs) associated with cardiac hypertrophy and fibrosis in the spontaneously hypertensive rat (SHR). In the current study, we tested the effects of Plzf gene targeting in the SHR using TALENs (transcription activator-like effector nucleases). SHR ova were microinjected with constructs pTAL438/439 coding for a sequence-specific endonuclease that binds to target sequence in the first coding exon of the Plzf gene. Out of 43 animals born after microinjection, we detected a single male founder. Sequence analysis revealed a deletion of G that resulted in frame shift mutation starting in codon 31 and causing a premature stop codon at position of amino acid 58. The Plzftm1Ipcv allele is semi-lethal since approximately 95% of newborn homozygous animals died perinatally. All homozygous animals exhibited manifestations of a caudal regression syndrome including tail anomalies and serious size reduction and deformities of long bones, and oligo- or polydactyly on the hindlimbs. The heterozygous animals only exhibited the tail anomalies. Impaired development of the urinary tract was also revealed: one homozygous and one heterozygous rat exhibited a vesico-ureteric reflux with enormous dilatation of ureters and renal pelvis. In the homozygote, this was combined with a hypoplastic kidney. These results provide evidence for the important role of Plzf gene during development of the caudal part of a body-column vertebrae, hindlimbs and urinary system in the rat

Translational regulation shapes the molecular landscape of complex disease phenotypes.

The extent of translational control of gene expression in mammalian tissues remains largely unknown. Here we perform genome-wide RNA sequencing and ribosome profiling in heart and liver tissues to investigate strain-specific translational regulation in the spontaneously hypertensive rat (SHR/Ola). For the most part, transcriptional variation is equally apparent at the translational level and there is limited evidence of translational buffering. Remarkably, we observe hundreds of strain-specific differences in translation, almost doubling the number of differentially expressed genes. The integration of genetic, transcriptional and translational data sets reveals distinct signatures in 3'UTR variation, RNA-binding protein motifs and miRNA expression associated with translational regulation of gene expression. We show that a large number of genes associated with heart and liver traits in human genome-wide association studies are primarily translationally regulated. Capturing interindividual differences in the translated genome will lead to new insights into the genes and regulatory pathways underlying disease phenotypes

Genetically Determined Folate Deficiency Is Associated With Abnormal Hepatic Folate Profiles in the Spontaneously Hypertensive Rat

Increased levels of plasma cysteine are associated with obesity and metabolic disturbances. Our recent genetic analyses in spontaneously hypertensive rats (SHR) revealed a mutated Folr1 (folate receptor 1) as the quantitative trait gene associated with diminished renal Folr1 expression, lower plasma folate levels, hypercysteinemia, hyperhomocysteinemia and metabolic disturbances. To further analyse the effects of the Folr1 gene expression on folate metabolism, we used mass spectrometry to quantify folate profiles in the plasma and liver of an SHR-1 congenic strain, with wild type Folr1 allele on the SHR genetic background, and compared them with the SHR strain. In the plasma, concentration of 5-methyltetrahydrofolate (5mTHF) was significantly higher in SHR-1 congenic rats compared to SHR (60±6 vs. 42±2 nmol/l, P<0.01) and 5mTHF monoglutamate was the predominant form in both strains (>99 % of total folate). In the liver, SHR-1 congenic rats showed a significantly increased level of 5mTHF and decreased concentrations of dihydrofolate (DHF), tetrahydrofolate (THF) and formyl-THF when compared to the SHR strain. We also analysed the extent of folate glutamylation in the liver. Compared with the SHR strain, congenic wild-type Folr1 rats had significantly higher levels of 5mTHF monoglutamate. On the other hand, 5mTHF penta- and hexaglutamates were significantly higher in SHR when compared to SHR-1 rats. This inverse relationship of rat hepatic folate polyglutamate chain length and folate sufficiency was also true for other folate species. These results strongly indicate that the whole body homeostasis of folates is substantially impaired in SHR rats compared to the SHR-1 congenic strain and might be contributing to the associated metabolic disturbances observed in our previous studies

Cap analysis of gene expression reveals alternative promoter usage in a rat model of hypertension

The role of alternative promoter usage in tissue-specific gene expression has been well established; however, its role in complex diseases is poorly understood. We performed cap analysis of gene expression (CAGE) sequencing from the left ventricle of a rat model of hypertension, the spontaneously hypertensive rat (SHR), and a normotensive strain, Brown Norway to understand the role of alternative promoter usage in complex disease. We identified 26,560 CAGE-defined transcription start sites in the rat left ventricle, including 1,970 novel cardiac transcription start sites. We identified 28 genes with alternative promoter usage between SHR and Brown Norway, which could lead to protein isoforms differing at the amino terminus between two strains and 475 promoter switching events altering the length of the 5' UTR. We found that the shift in Insr promoter usage was significantly associated with insulin levels and blood pressure within a panel of HXB/BXH recombinant inbred rat strains, suggesting that hyperinsulinemia due to insulin resistance might lead to hypertension in SHR. Our study provides a preliminary evidence of alternative promoter usage in complex diseases

Trans control of cardiac mRNA translation in a protein length-dependent fashion

Little is known about the impact of naturally occurring genetic variation on the rates with which proteins are synthesized by ribosomes. Here, we investigate how genetic influences on mRNA translational efficiency are associated with complex disease phenotypes using a panel of rat recombinant inbred lines. We identify a locus for cardiac hypertrophy that is associated with a translatome-wide and protein length-dependent shift in translational efficiency. This master regulator primarily affects the translation of very short and very long protein-coding sequences, altering the physiological stoichiometric translation rates of sarcomere proteins. Mechanistic dissection of this locus points to altered ribosome assembly, characterized by accumulation of polysome half-mers, changed ribosomal configurations and misregulation of the small nucleolar RNA SNORA48. We postulate that this locus enhances a pre-existing negative correlation between protein length and translation initiation in diseased hearts. Our work shows that a single genomic locus can trigger a complex, translation-driven molecular mechanism that contributes to phenotypic variability between individuals